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ABSTRACT 

Radio metric Jeep space navigation relies nearly 
exclusively upon coherent, i.e. 2-way, doppler and 
ranging for ail precise applications. These data 
t/pes and the navigational accuracies they can pro- 
duce are reviewed. The deployment of Hydrogen maser 
frequency standards and the development of Very Long 
Baseline Interferometry (VLBI) systems within the 
Deep Space Network is making possible the development 
of non-coherent , 1-way data forms that promise much 
greater inherent navigational accuracy. These data, 
closely paralleling the observables taken with VLBI 
arc much more sensitive to clock synchronizations, 
both time and frequency, and to instability during 
the measurement period itself than are the coherent 
data. The underlying structure between each data 
class and clock performance is charted. VLBI obser- 
vations of the natural radio sources are the planned 
instrument for the synchronization task. This method 
and a navigational scheme using differential measurements 
between the spacecraft and nearby quasars are described. 


I. INTRODUCTION AND BACKGROUND 

Navigation for deep space probes has always required the acquisition of 
long arcs of precise doppler data to determine accurately and reliably 
the spacecraft orbit prior to planetary encounter. The Doppler data 
used is taken in the so-called coherent, or 2-way, mode. In this 2-way 
mode, a ground station within the Deep Space Network transmits a stable 
frequency reference to the spacecraft which transponds a coherent ver- 
sion of the received signal back to the ground station for Doppler de- 
tection (Renzetti, ef al 1 and Melbourne^). At or near planetary 
encounter, the gravitational effects of the target body on the motion 
of the spacecraft are sufficient to produce a unique velocity profile 
in the Doppler that can be the dominant effect in determining where the 
spacecraft is. These effects usually arise too late, however, to be of 
benefit in the navigation process of determining and adjusting the 
planet-relative encounter conditions. The tracking which determines 
the planet-relative navigation accuracy is performed long before the 
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actual encounter and when the gravity forces on the probe are of little 
benefit in supplying the desired location information. Here the dominant 
information is supplied by the diurnal motion of the tracking stations 
as illustrated in Figure 1. A single pass of doppler tracking data will 
yield information about three quantities: 

1) The probe’s geocentric velocity by calculating the doppler 
shift with the diurnal effects removed. 

2) The right ascension, by determining the time of the meridian 
crossing via observing the null in the diurnal harmonic modu- 
lation. 

3) The cosine of the declination, by determining the amplitude of 
the harmonic modulation. 


The complete spacecraft state is determined by taking additional passes 
of doppler data on succeeding days and combining the information from 
each into a full solution. When available, ranging data greatly aids 
this process of determining the 6-dimensional spacecraft state to the 
requisite accuracy. In summary: 

1) range and doppler data give direct measurements of the line of 
sight probe distance and its time derivative, 

2) the diurnal modulation of the doppler yields an indirect 
measurement of the right ascension, a, and declination, 6, 

3) successive passes detect changes in a and <5 permitting indirect 
measurements of a, <5 and provide the data volume base for 
noise averaging and data consistency. 


Several detailed analyses of this situation have been performed - among 
-hem are: Hamilton and Melbourne , Curkendall and McReynolds 4 - giving 

the accuracy of the observable parameters as a function of doppler 
measurement accuracy, tracking system calibrations, probe geometry, 
length of the tracking pass, tracking station location errors, and ran- 
dom non-gravitational forces affecting the motion of the spacecraft. 


The accuracy to which the indirectly measured parameters can be estima- 
ted determine the overall accuracy of the complete orbit. Shown in 
Figure 2 is a plot of the system accuracy at several time points from 
the inception of the planetary exploration era to the present time - 
although the list is far from complete, the margin contains a tabulation 
of some of the technology improvements which enabled the accuracy evolu- 
tion as shown. The accuracy achievable from a given doppler and asso- 
ciated data calibre lion and processing system is seen to be a strong 
function of the nominal probe declination. This is because the estimate 
of declination is determined K observing the amplitude of the diurnal 
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harmonic signature, wr cos 5, and this amplitude becomes a stationary 
at 6 = 0. s 

The deterioration of the orbit performance at the lower range of declin- 
ation can be a serious inconvenience for inner planet exploration. For 
example, the allowable Viking Mars arrival date space was truncated due 
to declinations lower than 5 deg for arrival dates past September 1976. 
The performance versus declination shown in the figure thus manifests 
itself as a mission planning constraint. Fortunately, the tempo of the 
geometry evolution within the inner planet system is rapid enough so 
chat the periods of low declination and resulting loss of navigation 
ability are relatively short and can be usually tolerated. 

The situation can be much worse for outer planet exploration. For 
example, the Voyager I and II Saturn encounters are both such that 
neither spacecraft is above 5 deg declination during the 4 months prior 
to the critical planetary encounter. To meet this and like situations, 
the dual-station planetary ranging system has been developed. Unlike 
the doppler system, two tracking stations are involved as is shown in 
Figure 3. In order to measure the troublesome declination variable 
these stations need to be widely separated in latitude as shown. The 
difference of the two range measurements is proportional to B sin 6, an 
observable with the right structural relationship with declination; i.e., 
the measurement sensitivity to declination maximizes near zero rather 
than becoming stationary. The declination accuracy achievable can thus 
approximately be given by 


Bcos6 

where 0 A , is the overall accuracy of the system's ability to 
measure“£he range difference, Ap, and B is the polar projection 
of the baseline. Present mechanizations require that the two range 
measurements be performed sequentially, rather than simultaneously, with 
each station obtaining a measure of the round-trip light-time between 
the probe and itself. The measurements can then be referred to a common 
epoch by either modeling the probe's motion or accumulating doppler data 
during the intervening period. The accuracy of this system is currently 
limited to a performance of approximately as discussed by 

Christensen and Siegel^. This measurement accuracy, working with the 
baseline of the tracking stations at Goldstone, CA, and Canberra, 
Australia, yields a declination accuracy of approximately 1 ytad over 
the full normal operating range of declination. It has been an itera- 
tive process, but this performance and the current Voyager mission 
requirements are commensurate as discussed by Jordan. Figure 4 summar- 
izes the current performance of the radio metric tracking systems for 
both the doppler and the ranging measurements. As shown, the doppler 
system provides a performance at the 0.25 yrad level for high declina- 
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tions. This degrades slowly until 1 prad is reached; the dual station 
ranging system then provides a level 1 prad performance for the remainder 
of the declination space. 

II. VERY LONG BASELINE INTERFEROMETRY 


Since about the beginning of the current decade. Very Long Baseline Inter- 
ferometry (VLBI) systems have been developed in parallel with the coher- 
ent doppler and ranging systems just described. References 7-18 
which span this time period, describe the development and provide analyses 
necessary for a detailed understanding. In this paper we shall be con- 
tent to describe VLBI only in tutorial terms, emphasizing the similari- 
ties and contrasts between the VLBI data and the coherent range and 
doppler already discussed. 


In a typical VLBI system, each of two widely spaced antennas observes a 
single (broad band) radio source, e.g., a quasar, simultaneously record- 
ing the received signal over a specified frequency interval. The record- 
ings are digital in which the received voltage is digitized at the one 
bit level; timing information is added so that the recordings may be 
cross correlated later when brought to a central site. The (expected) 
cross correlation function can easily be shown to be approximately 
(Thomas-^-) : 


E f R(t , At) 



sin 77 WA t 
rr WAt 


cos 
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Where 


At = t - x 

g m 

t (t) = geometric delay as shown In Figure 5 

g 

X = a priori delay estimate inserted to bring the 

correlation function to near its maximum during data 
processing 

W = bandwidth of recorded signal 
4> i = At 

u>. = frequency at the center of the bandpass 

Using typical values of 2300 MHz and 2 MHz as the rf and recorded band- 
width respectively, this correlation function goes through one complete 
cycle for every change in At equal to the period of the rf frequency 
(<0.5 nsec). In addition it also manifests a sin x/x characteristic 
envelope reaching Its firsc null at 0.5 usee delay. These two compo- 
nents are called respectively "fast fringes" and the "delay function" 
(see Figure 6). For these same typical values, T (or equivalently .Ap/c, 
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Figure 5) can be measured directly, by adjusting T so as to maximize the 
delay function, to on the order of 10 nsec precision (3m in light-sec) . 
More powerful measurement;, can be obtained, however, in each of two con- 
ceptually different ways: 

1) Observation of the source continuously over the common 
visibility period of the two stations (approximately 4 hours on 
the baselines afforded by the DSN for sources near the ecliptic), 
produces a conti n ous record of the phase of the fast fringes 
versus time. I u record thus obtained will contain a diurnal 
sinusoidal modu ition term due to the Earth's rotation whose 
phase and av plit .de are parametric in the source location and 
baseline paranetcrs. This is exactly analogous to the single 
station coherent doppler tracking except that the equatorial 
baseline projection and longitude play the roles of the distance 
off the sp:..n axis, r , and station longitude, 9, respectively. 

The differential frequency of the two clocks replaces the 
geocentric veloc ty term observed by the coherent data. 

2) Observati .i of the source at a second center frequency, > 
produces a ;eccm< measurement of the fast fringe phase, 4> 2 > at 
a single instant of time. Then because 


9^ T g u> 2 - to 


a direct measurement of T can be obtained in the short time 
required to achieve a higM S/N for the (j> measurements (typi- 
cally 10 rain.) . This is the "bandwidth synthesis" technique, 
so called beca lse large effective bandwidths can be obtained 
without the ne d for commensurate high recording rates, and is 
widely used throughout the Vhol community (Rogers 9) . The 
geometry is exactly as for uhe differenced range measurement 
already discussed (cf. Figures 3 and 4); the T measurement 
obtained can be usee directly to estimate eithlr the baseline 
projection or the source location. With spanned bandwidths, 
l °2 ~ ^1’ otl t * ie order of 40 MHz, the precision of the measure- 
ment can easily be brought to the cm level; its accuracy is 
dominated by other effects such as clock performance and 
systematic calibration error'.. 

When estimating source lc 'tions with method two, a second 
baseline is usually cm{ iyed for the second component of 
position. An effective combination for the two baselines is 
to have a large do 1 >r component associated with the first and a 
large equatorial projection associated with the second so that 
they can produce largely uncorrelated estimates of a and 6. 
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Thus these two methods, often referred to as narrow-and wide- band VLBI 
respectively, have a one-to-one correspondence with the two coherent- 
modes, doppler and differenced ranging, normally used for spacecraft 
tracking. Their normal applications are duals of each other in that 
VLBI is usually employed to estimate the station baselines; precise 
source coordinates are needed to enable this. Coherent tracking is 
normally used to estimate the spacecraft coordinates; precise station 
locations are needed for this task. 

Although the VLBI and the coherent tracking modes each produce observ- 
ables with identical information content as just discussed, natural- 
source VLBI enjoys several inherent accuracy advantages over its coher- 
ent counterpart. These were discussed in some detail in a previous 
paper (Curkendall^) , but briefly they include: 1) wider bandwidth, 

2) more complete calibration of charged particles, 3) a ready means for 
calibrating electrical path delay variations in the station electronics, 
4) lack of significant proper motion in the natural sources themselves, 
and 5) freedom from needing to model the line-of-sight motion as is 
required in single station doppler tracking. 

There is a single major exception to the general advantages of the non- 
coherent data types - they suffer from a greater sensitivity to instabi- 
lity of the station master oscillator. This sensitivity and the compar- 
ison of it with that of the coherent data forms is treated in detail in 
the following section. 

III. MEASUREMENT ACCURACY VS. FREQUENCY STANDARD PERFORMANCE 

The two measurement classes just discussed, coherent measurements and 
non-coherent VLBI measurements differ markedly by the manner in which 
the station's master frequency standard departures from ideal enter and 
corrupt the measurements. 

In this section the four data types: 


\ Data 

^\£lass 
Data 

Type \ 

Two-Way 

Coherent 

Measurements 

One-Way 

Non-Coherent 

Measurements 

Narrowband 

Doppler 

(1) 

Narrowband 

VLBI 

(2) 

Wideband 

Differenced 

Wideband 


Ranging 

VLBI 


(3) 

(4) 


will each be analyzed and their sensitivity to clock performance charted. 
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Case 1 - Coherent Narrowband (Doppler) Data 


Consider the (highly) schematic diagram of a typical coherent doppler 
and ranging system implementation as shown in Figure 7. Counted, or 
integrated, doppler is obtained by broadcasting a stable reference to 
the spacecraft which coherently transponds the received carrier back to 
the same station for comparison with the original transmitted frequency; 
the difference or doppler frequency is integrated by means of a counter 
as shown. Assume for the purposes of illustration, that a unit step in 
frequency error occurs for a short period of time as shown in Figure 8 . 
This will enter the doppler extractor and be integrated to yield immed- 
iately a buildup of range error, Ap £ . It will also be transmitted to 
the spacecraft and return a round-trip light time, T, later and re-enter 
the doppler extractor, this time in the opposite sense - Ap £ will return 
to zero. If the doppler system has been tracking the spacecraft for T 
seconds, the accumulated effect of the time history of Af(t) is readily 
seen to be 


Ap_ 


c 

2f 


f 


Af (t) - Af(t-x) dt = 


r t 


2f L 

*'T-t 


Af (t)dt 


-f 

* / -T 


(T > T) 


Af (t)dt 


( 1 ) 


where c is the speed of light. 

That is, in a data point measured at T, frequency standard performance 
during the first and last x seconds counts, everything else cancels out. 
The factor of 2 appears so that Ap £ is .'ne error in the one-way range 
as measured by the two-way instrument. It is useful to design expres- 
sions which permit calculating Ap^ assuming Af is 
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i) a white noise process 

ii) linear with time 


i) White Frequency Noise 


With the white ncise assumption, the two integrals in (1) are clearly 
independent, the variance of Ap £ is then 



c 

2 



<5(u-v)dudv = 


j J2~h~~z 
2 v o 


( 2 ) 


where h 

o 


is the spectral density of the white frequency Af/f process. 


ii) Linear Frequency Drift 


Assume 


Af(t) 

f 


kt 


Then from (1), 0. = | k T T 

P E 


(3) 


where k is re-interpreted to be the standard deviation of the drift. 

Thus for frequency variations rapid relative to T, the error builds with 
/f and is independent of the tracking time; for frequency variations 
slow relative to T, the error builds as the product, Tx. 


Case 2 - Non-Coherent , Narrowband (VLBI) Data 


The non-coherent case is even more straightforward. Here the error in 
the measurement is easily seen to be proportional to the difference in 
the oscillators' frequencies at the two Earth-based stations integrated 
over the observational interval: 



( 4 ) 


where f = i ..tan' neous frequency at the ith station. 

The Af notation used in (1) is dropped here to emphasize that the mea- 
surement is sensitive to more than just the change in frequency over T 
or even T, it is sensitive to the "knowability" of the frequency diff- 
erence. The initial frequency offset is large enough so that in any 
VLBI experiment it must be considered an unknown and solved for - 
indeed, the determination of f^ (0) - f 2 ( 0 ) is often the reason for the 
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VLBI experiment. It becomes convenient, then, to re-define the problem 
and focus on a modified error term as 






Ap e 


i 

/ 


Af i - Af 2 


dt 


(5) 


where Af is understood to be the departure in the frequency from f (0) . 
Then 1 



where j is the familiar two-sample 


( 6 ) 


Allan Variance over smoothing time, T. Equation (4) is already in what 
is usually the most convenient form. For comparative purposes, however, 
it is interesting to recompute the effects arising from the white noise 
and linear drift models used earlier. 


i) 


White frequency noise 




6(u-v)dudv = 



(c.f.2) 


ii) Linear Frequency Drift 


(7) 


Af i - 


Af, 


~\]Z~ kt 
h k T 2 


Ap ' 


(c.f.3) 


( 8 ) 


Case 3 - Differenced Coherent (Ranging Data) 

A coherent range measurement is essentially a measurement of the round- 
trip light time itself , The clock error introduced in such a measure- 
ment is thus the absolute frequency error integrated over the light 
time. A differenced ranging measurement is sensitive to the integrated 
frequency differences, i.e.. 


Ap e “ 2f f f l “ f 2 du (9 ' 

J 0 

For precision differenced ranging measurements, this frequency differ- 
ence (or more precisely, the error in the knowledge of the difference) 
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must be held to within strict limits. For example, at the distance of 
Saturn, r % 10 1 * sec., must be < 3 xlO - ^ 3 for a Ap £ of less than 

.5 m. The implied frequency synchronization must be accomplished with 
traveling clocks or, as in more frequently the case, with a VLBI obser- 
vation session whose object is the determinations of 

The error, then, in applying the synchronization to a differenced range 
measurement can be roughly predicted as 



( 10 ) 


where At is the time between synchronization and application. This 
expression assumes that the synchronization operation has an associated 
error much less than (Af/f)^. This is not always the case and the VLBI 
method for synchronization brings up an interesting interplay between 
clock performance and its measurement. In an ideal experiment where all 
parameters influencing the interferometer phase (geometry, media trans- 
mission effects, etc.) are known save the clock offset, f - f_, itself, 
the clock performance during the experiment of duration, T, will contri- 
bute an error 


a 


f 


1 




( 11 ) 


where a. is given by (6) . 

Ap e 

Case 4 - Non-Coherent Wideband (VLBI) Data 


Here the measurement error is proportional to the clock time offset at 
the time of the measurement. Once again, this parameter must be mea- 
sured periodically in order that the knowledge of it can be held to 
within reasonable limits. The expression for the error using the same 
nomenclature as the first three cases would be 



( 12 ) 


where is understood to be the time of the last clock epoch synchro- 

nization operation. As before, if the synchronization operation is 
assumed accurate, the predicted standard deviation of (12) is readily 
seen to be 



where At = T - 
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Evaluation of (13) for even Hydrogen maser stabilities (Af/f 10 14 ) 
discloses that At cannot exceed 1/3 day if decimeter level measurements 
are sought. Because of this high sensitivity, it is common practice to 
include provision for both clock epoch and clock frequency synchroniza- 
tion integral to any natural source VLBI experiment. Equation (13) has 
relevance strictly only when envisioning a series of natural source 
measurements for clock synchronization whose results are applied to sub- 
sequent (or earlier) spacecraft tracking. 

Table I summarizes the four cases just discussed and repeats the four 
basic sensitivity equations. Table II is an attempt to tabulate the 
expected metric errors arising in the same four cases in terms of the 
two-sample Allan Variance. There are several approximations used in 
writing the expressions shown and this table should be viewed more as an 
ordered collection of the principles discussed here rather than a set of 
rigorous relationships. 

The non-coherent VLBI measurements are thus much more demanding of 
clock performance than their coherent counterparts. Indeed, the 
relative immunity of coherent doppler and range to clock variations en- 
abled reasonable performance at small light-time distances (<_ lunar 
distance) even with the crystal oscillators that were employed in the 
early 60's. With current rubidium standards (Af/f < 10 2 /day) good 
performance can be achieved throughout the terrestrial planet space. In 
contrast, VLBI errors are dominated by clock effects when rubidium and 
even cesium (Af/f < 10 13 /day) standards are used. The introduction of 
the Hydrogen maser with drifts better than 10~ 14 /day has in large part 
prompted the current interest in VLBI systems. At this clock performance 
level, the VLBI accuracy estimates given in the next section can be 
achieved and the non-coherent measurement class can successfully compete 
with its coherent counterpart. 


IV. VLBI AS A NAVIGATION TOOL 

The maturing VLBI technology can be applied to the spacecraft navigation 
problem in each of two conceptual ways: 1) calibration of the DSN for 

use as an otherwise conventional radiometric network, and 2) direct 
spacecraft signal tracking with the VLBI data acquisition and processing 
systems . 

VLBI for DSN Calibrations 


In July of '79, the first operational VLBI system (as contrasted with the 
R&D systems which have produced the results discussed above) will begin 
taking routine measurements on the California-Spain and California- 
Australxa baselines. Initially this system, operating in conjunction 
xtfith a developed precision source catalog, will have the capability of 
operating with an overall accuracy at the 0.05 urad level. This capa- 
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bility will be exploited to calibrate the UT and PM variations, refine 
the knowledge of baseline vectors, and determine the relative station 
clock frequency offsets. The UT/PM calibration improves the knowledge 
of the effective station location coordinates for both range and doppl 
the baseline solutions directly improve the differenced range data. 

Both of these error sources will be controlled to the 50 cm level as 
compared to their current 1-2 m effective levels. The interstation 
clock frequency calibration will be accurate to better than 1 x 10 13 , 
bringing the clock's contribution U *'he differenced range error 
(Eq. (10)) to well under 50 cm even at 10 AU. 

Direct VLBI Navigation 

The spacecraft can be treated as another VLBI radio source, albeit one 
with proper motion. Present spacecraft do not emit signals with band- 
widths wide enough to really be considered wideband sources, however. 
Instead, the power is centered within a few MHz of the carrier, effec- 
tively enabling tracking at only the carrier center frequency so that 
the fast fringes may be observed but bandwidth synthesis at beneficial 
accuracies cannot be performed. Thus, only narrowband VLBI can be made 
available for already launched spacecraft. Save the one feature of the 
wide bandwidth, all of the inherent advantages of the non-coherent 
data class can be capitalized on, however. 

The favored operational procedure in narrowband spacecraft VLBI is to 
track the spacecraft in close conjunction with the tracking of a nearby 
quasar ('v <5 deg away in angular measure) whose source location has 
already carefully been determined. In this procedure, known as AVLBI, 
the two antennas track the spacecraft for a few minutes and then quickly 
align on the quasar. This is repeated for the entire duration of the 
'< hour tracking pass and the fringe phase history differential between 
the spacecraft and the source is constructed. The differences in a and 
6 between the two sources are then estimated. This procedure has the 
advantage that most of the systematic error sources, being nearly common 
to both the spacecraft and quasar phase histories - baseline coordinates, 
instrumentation calibration, neutral and charged particle media effects, 
and clock imperfections - are diminished by the differencing operation. 
The precise natural source location estimate is quickly transferred to 
the spacecraft by this technique. 

It is relatively simple to alter the broadcast spacecraft waveform to 
make it suitable for wideband VLBI tracking. Modulating a high fre- 
quency subcarrier with a signal generated from a noise diode could 
produce noise channels suitable for bandwidth synthesis with character- 
istics nearly identical to that being received from the quasar. A far 
better technique is to transmit a specific ranging code generated on- 
board the spacecraft and c. rrelate the received spacecraft signal at 
each station against a locally generated model of that same code in a 
manner similar to that used in a conventional 2-way ranging machine. 
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In contrast to the current ranging technique, the signal would not be 
demodulated and tracked via the closed loop receiver, however. Rather, 
the spacecraft signal would follow the same rf chain as the quasar 
signal, phase calibration signals would be introduced, and detection 
would take place after digitization so as to preserve the near perfect 
commonality with the natural radio source tracking essential for careful 
system calibration. Each station executes this procedure, a one-way 
range is calculated, and the difference taken. This difference contains 
the differential station clock epoch error which must be carefully cali- 
brated with conventional VLBI. This calibration, along with UT/PM 
calibration, can be performed via either an earlier multiple-source VLBI 
run and applied to the Differential One-Way Ranging (DOR) data or the 
DOR data can be obtained in conjunction with the track of a single very 
nearby quasar in a manner similar to the narrowband AVLBI method already 
( described. Which method is employed operationally would be determined 

by the availability and strength of a nearby quasar, the number of 
spacecraft that need to be tracked, the availability of transmission 
lines for the VLBI data, the time criticality of the spacecraft 
results, and the ability of the clocks to hold an epoch synchronization. 

An extensive error analysis and discussion of the wide-band VLBI, or 
DOR, technique was presented by Melbourne and Curkendall in an earlier 
paper^O. As a result of that and subsequent analysis, it is felt that 
a confident prediction of the performance of such a system can be made 
at the 0.05 prad level. The narrowband AVLBI can perform at similar 
levels for spacecraft at high declination, but would degrade with the 
1/tan 6 characteristic inherent for all narrowband tracking. 

These accuracies, along with the accuracy of the coherent tracking 
earlier discussed, is summarized in Figure 4 . In contrasting the wide 
and narrowband approaches, the constant performance of the wide-band 
method with declination is of the utmost importance. In addition, the 
tracking time required for a complete observation is much smaller for 
the wide-band system since no diurnal signature need be observed. Ob- 
servations from two baselines are required however. The narrowband 
approach suffers from a greater sensitivity to systematic error sources 
so only Lne AVLBI mode should be used to meet a precision requirement. 
The wideband is more immune and offers greater tracking strategy flex- 
ibility. It is thus less dependent on the existence of a nearby quasar. 
This could be important for some applications since at the present 
time,, few suitable sources have been found in the portion of the 
celestial sphere where the galaxy partially obscures extragalactic 
observations (Preston et al^l) . This results in two bands along the 
ecliptic plane, each measuring about 30 deg in extent, nearly devoid of 
suitable sources. Finally, the post correlation data processing is 
more difficult for the narrowband AVLBI data. The key to the accuracy 
is the accurate construction of the accumulated phase delay change over 
the 4 hour period for both the spacecraft and the quasar. Care must be 
taken to ensure that the phase is extrapolated properly to within one 
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rf cycle (13 cm at S; 3.5 cm at X-band) during the time periods when the 
antennas are moving or are on the opposite source. 

On a more positive note, narrowband AVLBI can be very useful for track- 
ing existing spacecraft and offers a greater sensitivity for the detec- 
tion of proper motion over short time periods (< 4 h) . 


1 


Sensitivity of Differential VLB1 to Clock Performance 

Differential-, or A-VLBI, is intermediate in clock sensitivity between 
the coherent or non-coherent data forms discussed in Section III and 
for that reason deserves special treatment in this concluding sub-sec- 
tion. 


This sensitivity is generally proportional to the angular separation 
between the natural source and the spacecraft. For example, in narrow- 
band AVLBI, when the sources are both within a single beamwidth of the 
antennas, the normal switching between the two as described above is not 
necessary, both sources are tracked all the time and the effect of 
clock variations cancels out. All that is needed in this case are 
clocks with coherence times long enough to detect the signals ('v a few 
minutes). Beyond a single beamwidth separation, antenna switching is 
necessary and the individual spacecraft and quasar records will have 
complimentary data outages as shown in Figure 9. For separations just 
greater than the single beamwidth limitation the data differencing oper- 
ation would difference the two data streams as close together in time 
as possible or one full cycle time, denoted as At in the diagram. 
Consider for the moment, that the clock instability is the sole error 
source, t^e measured Ap for the quasar would be: 


Ap (T) = Ap (T) + 
mq q 


tf„ T £ i - 


dt 


(14) 


An identical expression for the spacecraft measurement, Ap , (T) can 
also be written. Their difference with the time shift. At is / 


AAp (T, At) = 

in 


It is useful to calculate the expected rms value of the clock induced 
error for the white frequency noise and linear drift examples used 
earlier. 


Ap (T + At) - Ap . (T) 
mq m s/c. 


Ap - Ap . + 

q s/c 


f t+t 


(15) 


f l ' f 2 dt 
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i) 


White Frequency Foise 


a.. = cj 2h At (c.f.2&7) 

AAp v o 


( 16 ) 


ii) Linear Frequency Drift 

°AAp = T (At2 + 2TAt) (c.f .3&6) (17) 

In comparing these expressions with those derived earlier, note that if 
the switching time is shorter than the spacecraft round trip light time 
(At can be on the order of 10 minutes, ~ is often several hours), the 
AVLBI data is actually less sensitive tu the white frequency noise than 
is two-way counted doppler. The expression for the linear drift model 
is a good illustration that the sensitivity to more systematic clock 
errors is intermediate to that of the coherent and non-coherent forms 
discussed in Section III. 

For larger source separations, the situation is somewhat more complex. 

As the angle increases, -t soon becomes apparent that if both sources 
are tracked over the same time periods (e.g. t = 0 to T) except for the 
alternating outages, the cancellation of other error sources will not be 
nearly so complete as if better strategies had been employed. For 
example, suppose that the sources have equal declinations but differ in 
their right ascensions and that the dominant error source is tne tropo- 
spheric scale height. It should be clear that an this circumstance, the 
best strategy would be to track each source over the same range of hour 
angles and stagger the time intervals as shown in Figure 10. 
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In this circumstance, the effective At would grow to a much larger value 
than is strictly needed for the switching time operation as is illuscra- 
ted in the figure. More generally, the global tracking scenario and the 
effective At must be chosen to minimize the sum of all the error sources, 
a problem beyond the scope of this short discussion. Once chosen how- 
ever, the clock’s contribution to the accumulated error can be written 
directly in terms of the Allan Variance as 



(18) 


In wideband AVLBI, the tracking strategy is simplified to a single obser- 
vation of each source. The clock's contribution to the error in the 
difference of the two resulting time delay measurements are 1) the error 
suffered internal to the measurement itself, and 2) the clock offset 
drift in between the two measurements, spaced At apart. The determina- 
tion of the first of these is beyond the scope of this discussion (but 
shou'd be small), the second is given by (18) just as in the narrowband 
cas°. The real difference between the narrow and wideband cases is that 
the narrowband transformation to right ascension and declination esti- 
mates is itself more sensitive to a time delay error change (by about a 
factor of 5 at high declination, growing gradually worse as declination 
is reduced) than is the wideband transformation sensitivity to time 
delay error. 

A single numerical example should serve to put these relationships in 
perspective. During the 1979 Voyager encounters with Jupiter, a narrow- 
band AVLBI demonstration is planned whose accuracy goal is set at the 
.05 ]jrad level. A single quasar, OJ287, is being used as the reference 
natural source throughout the several month experimentation period; the 
appropriate At is as large as 45 min. The error control needed to 
achieve this accuracy is approximately 7 cm of range change error during 
the 4 hour integration period. If 1 cm is the clock's allocation, 
the two-sample Allan Variance clock performance required is (from (18): 

If a.. = .01m 

mp € 

At - 45 min. 


Af 

f 


a AAp e 

A c At 


8.8 x 10 


-15 


If wideband AVLBI were possible, the same accuracy level could be 
achieved with a total error budget of about 40cm. If 1/7 of this were 
allocated to the clock, the Af/f specification could be relaxed to 
5 x 10 14 . Alternately, and probably of more practical importance, a 
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clock operatin' at 8.8 x 10 15 could be used to lengthen the time 
between epoch calibration and spacecraft use; i.e., the 45 min. 
could be stretched to several hours. 
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TABLE I: Comparison of Data Type Sensitivity to Station Master Clock Errors 
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: instantaneous frequency of i cn station clock 
; Tracking Pass Duration 
! Round Trip Light Time to S/C 

All errors are given in range difference or integrated range-rate (m) . 
















Table II. 

Approximate Evaluation of Expected Metric Error, a. , In 

/Af\ 2 Ap e 

Terms of the Two Sample Allan Variance l-j-l • 
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Variance 
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c = speed of light 

= t ' w ° sample Allen Variance evaluated at t smoothing time 

t = round trip light-cime to coherent transponder 

T = time from beginning of Ap integration 

^f = error from frequency synchronization operation 
s 

Ir = error from clock epoch synchronization operation 


At ~ time from clock synchronization to T 

© , an operator impling addition of errors in the rms sense 
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Figure 1. The Primary Information Source for Pre-Encounter Navigation Arises from the 
Diurnal Motion of the Tracking Station. 
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Figure 2. Doppler System Performance versus Time and Spacecraft Declination 
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Figure 3. The Difference of Two Ranging Measurements 
from Stations Widely Separated in Latitude 
can Accurately Measure Declination of Probe 
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Figure 5 
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Accuracy Performance of Doppler and Ranging System versus 
Nominal Probe Declination. 



Wideband VLBI can Estimate Source Location by Measuring 
the Differential Time of Arrival. 
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Figure 6. Band-limited White 

Noise VLBI Correlation FAST FRINGES 




Figure 7. Coherent Doppler, Range Cround Tracki ,ig System 
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Figure 8. Counted Doppler Response, Ap,_, To a Unit Step 
of Af/f. K ~ } 



Figure 9. Narrowband AVLBI Phase Record Showing Alternate 
Data Gaps Due to Antenna Switching. 
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Figure 10 


Narrowband AVLBI Tracking Scenario 
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Questions and Answers 


DR. TOM CLARK, NASA Goddard Space Flight Center: 

I might update a couple of comments that were made in that one. I 
believe there has been one successful spacecraft-to-quasar delta 
VLBI experiment which was done by Shapiro et al. at MIT, and Newhall 
et al. in your shop, on the differential position of the Viking 
Oroiter, and the OJ-287 quasar. That was done about a year ago, and 
the data is still being processed, but it appears that it was 
successful . 

I just wanted to stress one thing in terms of the angular mea- 
sures that were mentioned here. The units, milliseconds of arc and 
so forth, were stressed a few times, and I just wanted to remind you 
what a millisecond of arc was. 

If I took this quarter and gave it to Dave, and told him to 
take it back to Pasadena and tried to look at it from here, that is 
about a millisecond of arc. That is half a nanoradian. 

I would point out that there is at least one pair of quasars 
which have been done to a few tens of picoradians in terms of 
differential positions by a delta VLBI type technique. They are 
only a degree apart on the sky, but it is the size, roughly, of 
George's eyeball on this thing, in the same analogy as before. 

MR. CURKENDALL: 

Yes. Not only has there been one experiment with the Viking space- 
craft, there are something like 20, an unknown number of which are 
successful, sitting in our data hoppers. 

By my remark, I meant I don't think anybody has been ever able 
to run an experiment where the accuracy had to be seven centimeters, 
and then be able to prove that it was. And I think that is accu- 
rate; it has probably never been done. And I think this is the same 
thing whether you can really wring that kind of performance out of 
the system or not. 

DR. CARROLL ALLEY, University of Maryland: 

Unfortunately, I arrived late for your talk. You may have mentioned 
this at the beginning. But it is worth pointing out, I think, that 
with this coherent tracking, and with a good transponder on these 
deep spacecraft, that there is a potential of measuring low fre- 
quency gravitational waves that may well be occurring in the uni- 
verse. That is, the part in 10^ that you are striving for is about 
three orders of magnitude too insensitive, according to the current 
estimates of my friend Kip Thorne and o'her people. 

Nevertheless, I would like to submit, as in an earlier discus- 
sion today, that we do not know everthing about the universe, and 
that there may well be gravitational radiation of higher amplitudes 
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than is predicted, and I hope that you will be keeping a very care- 
ful watch, even at Lhe part in 10^ level. The strain induced, the 
delta L over L in the distance between the spacecraft and the earth, 
is the same as tne strength of the g^aviational radiation. 

MR. CURKENDALL: 

1 understand. 

DR. ALLEY: 

And so there may just be some relic radiation of this amplitude 
around that would show up in such measurements. 

MR. CURKENDALL: 

Yes. I think your word "relic" is the key here. The thing I don't 
like about the Kip Thorne calculations is it goes something like 
this: If you assume the collapse of something like 10^ solar 
masses, you will get a differential movement between the spacecraft 
and the earth. 

And maybe this occurs something like once every 30 years. You 
will get a differential change in length between the earth and the 
spacecraft of about one and a half mil.imeters; and the problem of 
sitting around for 20 years waiting for that to happen is immense. 

What seems more feasible to me is to look at what you have just 
said: Sit there and look for the background radiation. If you can 
put an X-band uplink on the spacecraft, and get masers down to a few 
parts in 10l 5 , the numbers already work out that you should be able 
to see enough gravitational energy that would close the universe. 
And that is much better than waiting around. 

And the advantage of that experiment is if you go out and you 
look and you don't see anything, you can go home. And that is 
really important. 

DR. VICTOR REINHARDT, NASA Goddard Space Flight Center: 

Just one question and it is a partial comment: Is there any possi- 
bility of using telemetry information to increase the bandwidth and 
narrowband VLBI to give you an effective wideband? 

MR. CURKENDALL: 

The telemetry, right now, runs about 3o0 kilohertz subcarriers and 
at around 100 kilobits a second. You see, relative to the earth's 
satellites, where the signal and the noise are so much better, the 
telemetry rate coming back is not that high. And that is not very 
much bandwidth spreading. 

We are doing that as part of this demonstration program. In 
fact, we have a mode where we leave the subcarrier on, but turn off 
the modulation so you get the nice pure sine waves, and you look at 
the harmonics of those sine waves. 





It is just like the Goddard side-tone ranging system, in terms 
of being able to detect that. 

But we detect it through the open loop RF chain, and calibrate 
with the quasar signals. And so, yes; you can do it as a restricted 
bandwidth sort of thing. 

What we hope to put on Galileo is tones of plus-minus 20 mega- 
hertz on down, and that ought to do it. 

MR. DAVID W. ALLAN, National Bureau of Standards: 

I would think that this level of sensitivity, you would be sensitive 
to lunar crustal tidal movements of the mantle of the earth. This 
was not mentioned. I was just wondering if this is a problem. Is 
this not of the order of a meter or two? I am not sure. 

MR. CUR KENDALL: 

Yes, they are; and I think our Chairman can give better words on 
that. I think the bottom line is that you have to model but it 
is thought to be model able. 

DR. TOM CLARK, NASA Goddard Space Flight Center: 

Yes, there are a number of these little subtle, insidious effects 
that affect the stations on the earth. The particular one with 
tides is not as bad as many of the others because the tides have a 
semi-diurnal signature, and it is very easy to pull semi-diurnal 
signatures out, because all of the things that Dave was talking 
about were diurnal signatures. It is those diurnal ones that are 
much more insidious and give you a lot more trouble; things like 
diurnal polar motion, for example. 

MR. CURKENDALL: 

Do you know what you are most sensitive to? I once did a spectrum 
analysis, and for reasons I have never quite understood, you are 
really most sensitive to twice diurnal rates. With a given amount 
of power, you are more sensitive to something with a 12-hour period 
than a 24-hour period. I have never quite understood that, physi- 
cally. 
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